Effect of phase fluctuations on INS and NMR experiments in the pseudo-gap regime 

of the underdoped cuprates 
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We present a theory for inelastic neutron scattering (INS) and nuclear magnetic resonance (NMR) 
experiments in the pseudo-gap regime of the underdoped high-T c cuprates. We show that su- 
perconducting phase fluctuations greatly affect the temperature and frequency dependence of the 
spin-susceptibility, x", probed by both experimental techniques. This result explains the appear- 
ance of a resonance peak, observed in INS experiments, below a temperature To > T c . In the same 
temperature regime, we find that the 63 Cu spin-lattice relaxation rate, 1/Ti, measured in NMR 
experiments, is suppressed. Our results are in qualitative agreement with the available experimental 
data. 



PACS numbers:74.25.-q,74.25.Ha,74.25.Jb,74.40.+k 

Over the last few years intensive research has focused 
on the origin of the pseudo-gap region in the under- 
doped high-T c cuprates ^Hffl- This part of the phase 
diagram, below a characteristic temperature T» > T c , 
is characterized by a suppression of the low-frequency 
quasi-particle spectral density, as observed by angle- 
resolved photo-emission (ARPES) || and scanning tun- 
neling spectroscopy (STS) experiments || . For the same 
compounds, inelastic neutron scattering (INS) experi- 
ments have revealed a sharp magnetic mode, the reso- 
nance peak, below T* , in contrast to the optimally doped 
cuprates, where it only appears below T c |Hj. Moreover, 
nuclear magnetic resonance (NMR) experiments find a 
strong decrease of the 63 Cu spin-lattice relaxation rate, 
1/Ti, below T, 

These experimental observations put tight restrictions 
on the proposed theoretical scenarios for the pseudogap 
ascribing it to spin-charge separation JtJ, SO (5) symme- 
try H , condensation of performed pairs J5 [_ and spin- 
fluctuations ||]. Emery and Kivelson (EK) [jj proposed 
that, due to the small superfluid density of the under- 
doped cuprates, thermal fluctuations in the phase of the 
superconducting (SC) order parameter destroy the long- 
range phase coherence in the pseudo-gap regime, while 
preserving a finite local amplitude of the order param- 
eter, |A(r)|. In this communication, we argue that the 
presence of phase fluctuations provides an explanation 
for the results of INS and NMR experiments discussed 
above. We show that these fluctuations greatly affect 
the temperature and frequency dependence of the spin- 
susceptibility, x", probed by both experimental tech- 
niques. Support for the existence of phase fluctuations 
comes from recent high frequency transport experiments 
by Corson et al. |TJ]. They demonstrated that the SC 
transition in underdoped B^Si^CaC^Os-a (Bi-2212) is 
of the Kosterlitz-Thouless (KT) type, where at T c = 7AK 
the unbinding of thermally excited vortex-anti-vortex 
pairs destroys the long-range phase coherence. However, 
they also concluded that |A(r)| vanishes at a tempera- 
ture Tq ~ 100K, while the onset temperature, T», for 



the pseudo-gap regime is much higher. For this reason 
we will focus our analysis on the region T c < T < T$. 

The starting point for our calculations is the mean-field 
BCS Hamiltonian in which the phase, 0(r), of the super- 
conducting order parameter, A(r) = \A(r)\ e tS ^ r \ varies 
on the scale of the phase coherence length, £g. Such a 
nonuniform phase 6(r) is treated via a gauge transfor- 
mation 



(1) 



where is the creation operator of the original electrons. 
This transformation induces a coupling of the ^-fcrmions 
to a local superfluid flow v s (r) = V0(r)/2m, (we set 
h = 1) whose thermodynamic properties are determined 
by the 2D-XY Hamiltonian 



Hxy = K (T) 
k B T 2 



J d 2 r|V6»(r) 



(2) 



where K (T) = n s (T) / (Amk B T) is the "bare" phase stiff- 
ness and n s (T) is the 2D superfluid density per Cu02 
layer which for a d-wave superconductor is given by 
n s (T) = n s (0)(l - T/T ), where T is the BCS mean 
field temperature 0. In order to compute \" m the 
presence of phase fluctuations we first compute it for a 
given configuration of v s (r) and subsequently perform a 
thermodynamic average over the ensemble specified by 
Eq.(||). Our approach is similar to the one adopted by 
Franz and Millis (FM) who computed the single parti- 
cle Green's function, G(k,uj), in the pseudo-gap regime. 
They showed that the quantity which determines the 
ensemble average is the correlator W = m 2 Vp(v 2 ) /2, 
whose temperature dependence they extracted from fits 
to ARPES and STS experiments. In the following we 
show that W(T) can also be obtained from the experi- 
ments by Corson et al. |l3| 

Assuming that the superfluid velocity is purely due to 
transverse phase fluctuations, we have 



W(T) = TT 2 V 2 F J 



rf 2 qG(q) 
47T 2 q 2 



(3) 
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where G(q) = (n q n_ q ) is the vortex density correlator. 
In the limit of large vortex density, this correlator is given 
by [0 G(q)- 1 - ^ 2 K {Q + q~ 2 ). Evaluation of the 
integral in Eq.(||) with wave- vector cutoff A = 27t£q L 
yields 
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Here, the phase coherence length is given by £^~ 2 (T) = 
ATT 2 Ko(T)np(T), where np(T) is the density of free vor- 
tices, and we used the BCS result £gl — vp/(ttA ). 
Corson et al. found in their analysis that np(T) — 
(2A/TT^ L )exp(-8CK (T)) with A and C being con- 
stants of order 0(1). Assuming that the functional form 
of rip(T) (and thus that of W(T)) remains the same for 
all underdoped Bi-2212 compounds, with A and C be- 
ing the only doping dependent parameters, we present 
W(T) from Eq.(|j) in the inset of Fig. [| together with 
FM's fits §] to STS experiments in Bi-2212 (T c = 83K) 
§. With A = 0.1 and C = 0.6, we find good quanti- 
tative agreement of our theoretical results with those of 
FM up to To « 150 K. At this temperature, the above 
approximation, as well as the analysis by FM, presum- 
ably break down since y/W(T) becomes of the order of 
the maximum superconducting gap. 

We now turn to the appearance of the resonance peak 
in the pseudogap region. Morr and Pines (MP) Jl5| re- 
cently argued that the resonance peak in the supercon- 
ducting state arises from a spin-wave mode whose dis- 
persion is given by 



= A 



2 

sw 



|q-QI 
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where A sw is the spin- wave gap, c sw is the spin- wave 
velocity and Q is the position of the magnetic peak in 
momentum space. Starting from a spin-fermion model 
fl6|| , MP showed that this mode is strongly damped in 
the normal state, but becomes only weakly damped in the 
superconducting state, if A sw is less than the gap, u) c , for 
particle- hole excitations with total momentum Q. These 
excitations connect points on the Fermi surface (FS) in 
the vicinity of (0, it) and (it, 0) ("hot spots"), and thus for 
a d-wave gap Ak = Ao (cos(fc x ) — cos(fc a ))/2, u c as 2Ao. 
MP computed x using the Dyson-equation 



X 



Xo 



n 



(6) 



where xo 1S the "bare" susceptibility and II is the bosonic 
self-energy given by the irreducible particle-hole bubble. 
For xo, MP made the experimentally motivated ansatz 



Xo 



(7) 



where u> q is given in Eq.(|^). In the superconducting state, 
one obtains for II to lowest order in the spin-fermion 
coupling g 



II(q,*w n ) = ~9 2 TJ2 {G(k,in m )G(k + q,in 

k , 771 

+F(k, itt m )F{k + q, itt m + iuj n ) J , (8) 

with G and F being the normal and anomalous Green's 
functions. Since, within the spin-fermion model, xo is 
obtained by integrating out the high-energy fermionic de- 
grees of freedom, it is largely unaffected by the onset of 
superconductivity or the pseudo-gap. Moreover, MP ar- 
gued that due to fermionic self-energy corrections, Re II 
in the SC state only leads to an irrelevant renormaliza- 
tion of A sw and c sw . Since the same argument also holds 
within our scenario for the pseudo-gap region, we neglect 
Re n in the following. 

On the other hand, Imll which determines the damp- 
ing of the spin excitations, changes dramatically in the 
SC state due to the opening of a gap in the fermionic dis- 
persion. Consequently, we expect phase fluctuations to 
strongly affect Imll. Moreover, in each polarization bub- 
ble present in the RPA expansion of Eq. @ , the electron- 
hole pairs probe different parts of the sample and thus 
independent configurations of thermally excited super- 
currents. It then follows that the susceptibility, Xp/j m 
the presence of phase fluctuations is obtained from Eq. (||) 
by using lmll p f averaged over the thermodynamic en- 
semble determined by Eq.(|2|). 

Before we discuss the effect of phase fluctuations on 
Imll, we shortly review its form in the normal and SC 
state. Extending Eq.(@) to the normal state, we obtain 
Imn(Q) = Ag 2 Lo / {-KVp) jl7|], where vf is the Fermi ve- 
locity at the hot spots. In contrast, in the SC state, in 
the limit of T oj c , we find to order 0(T/cu c ) 



ImII(Q,w) = 



Ag 2 uo c 



E{^/l-uP)6{oj-l) 



TTVf 

^1^ + 1)6^-1) , 



(9) 



where 6{x) is the Heavyside step function, E{x) is the 
complete Elliptic integral of the first kind and Cu — u)/u> c . 
Thus, Imll vanishes for frequencies below u) c . In Fig. |l| 
we present the frequency dependence of Im n in the nor- 
mal and SC state. 

We now consider the effect of phase fluctuations on 
Imll. Note that n, Eq.(^), and thus x, Eq.dj|), are in- 
variant under the gauge transformation, Eq.pl), in con- 
trast to G(k, u>), considered by FM. Thus G,F can be 
straightforwardly calculated using the ^>-fermions. In the 
limit kp^g 3> 1, where kp is the Fermi momentum at the 
hot-spots, the interaction of the ^-fermions with the su- 
perfluid flow leads to a Doppler shift in the ^-excitation 
spectrum B,H| given by 



±£>k 



(10) 



where £k is the fermionic dispersion in the normal state, 
Dk = mvpfk) • v s is the induced Doppler-shift and Vp(k) 
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FIG. 1. Imll in the normal (black line) and SC state 
(dashed line), and IniIIp/(a;) in the pseudo-gap region for (a) 
W = 0.015 (dotted line), and (b) W = 0.1 (dashed-dotted 
line) . 



is the Fermi velocity. In the limit T <C -Dk "C Aq, Imll 
for a given superfluid velocity is obtained from Eq. via 
the frequency shift 



{D X + Dy 



(11) 



Similar to the case of the fermionic spectral function , 
the thermodynamic average of Imll over the ensemble 
specified by Eq.(|^) is obtained by convoluting Imll with 
a Gaussian distribution of Doppler shifts of the form 



P(D a ) = 



1 



\/2nW 



exp 



D 2 
2W 



(12) 



where a = x,y. In the limit yW <C T <C uj c , we can 
perform this convolution analytically and obtain 



Im H p f(io) 



2v% 



(l + O) 



W(T) /tt exp -i= 




where $(x) is the error function. It follows from Fig. [I], 
in which we present the spin-damping for two different 
values of W — W/Aq, that the effect of phase fluctuations 
on IrnlT is two-fold. First, they lead to a non-zero value 
of ImHpj(u;) for to < lu c , in contrast to the form of Imll 
in the superconducting state where the spin-damping at 
T = vanishes below uj c . Second, the spin-damping 
below lo c increases with increasing W while at the same 
time, the sharp step in Imll is smoothed out. Note that 
in the pseudo-gap region, r«w c , and consequently the 
temperature dependence of lmH p f(uj) is determined by 
that of W(T). 

Finally, inserting ImILpj(a;) into Eq.(^|), we obtain 
Xp/(Qi w ) m the pseudo-gap region. In Fig. || we present 
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FIG. 2. The resonance peak in the pseudo-gap region for 
S^f- = l/(70meV),A sro = 35meV, lj c = 65meV and W(T) 

from the inset. Inset: W(T) from Ref. (points) and from 
EqM with A=0.1 and C=0.6 (solid line). 



our theoretical results for the frequency and temperature 
dependence of the resonance peak in Bi-2212 (T c = 83K ), 
using the W(T) shown on the inset. We find that, as the 
temperature is increased above T c , the resonance peak 
becomes broader, while its intensity diminishes. Since 
W(T) is a monotonically increasing function of tem- 
perature, it follows from Fig. |l| that the spin damping 
for lu rj A sw < u c also increases with temperature, 
giving rise to the behavior of the peak intensity /width 
shown in Fig. ^|. Unfortunately, no experimental data 
for the temperature dependence of the resonance peaks 
in the pseudo-gap region of underdoped Bi-2212 are cur- 
rently available. However, our results are in qualitative 
agreement with the experimental data on underdoped 
YBa 2 Cu 3 6+:!; @. 

We now turn to the second experimental probe of x", 
the 63 Cu spin-lattice relaxation rate, \/T\. For an ap- 
plied field parallel to the c— axis, 1/Ti is given by 



1 
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where 



F c (q) = [A ab + 2B (cos(q x ) 



(14) 



(15) 



and A a b and B are the on-site and transferred hyper- 
fine coupling constants, respectively. The spin-lattice re- 
laxation rate in the mixed state, i.e., in the presence of 
a superflow, was recently considered by Morr and Wor- 
tis (MW) [p| . Using the low- frequency limit of Eqs.(||) 
and (||), they found that the temperature dependence of 
1/Xi is determined by the set {D n /T}, where D n is the 
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FIG. 3. 1/TiT in the pseudo-gap region of underdoped 
Bi-2212 (Tc = 79K). Solid line: theoretical fits with 
A=0.05 and C=0.3. Filled squares: experimental data taken 
from Ref. p2|. We assumed a constant background factor 
- - 0AK~ V^ 1 



Doppler-shift at the nth node (see Eq.(lO)). In the limit, 
\D n /T\ > 1, they obtained 



— = -V^(q 



C 

N 



(16) 



where C = (k B /^(ag^je) 2 / (4v F v A ) 2 , v A = |<9A k /<9k| 
at the nodes, and 



Q 



2)2 



(17) 



Here, q^- is the wave- vector connecting the nodes i and 
j, and £ is the magnetic correlation l eng th. In the limit 
T < yJW{T), the convolution of Eq.@ with the Gaus- 
sian distribution of Eq.(jl^) can be performed analyti- 
cally, and we obtain 



T X T J 



= /3W(T) 



(18) 



Pf 



where (3 = 4C(.F(0) + -F(qi, 3 ) + £ The con- 
stant P can be experimentally obtained |19| by fitting 
(TxT)- 1 at T < T c with the d-wave BCS expression 

(TyT)- 1 = (3^-T 2 . Note that the relaxation rate in 
Eq. ( |l8| ) directly reflects the strength of the classical phase 
fluctuations. In Fig. || we present our theoretical re- 
sults for {T x T)~j, Eq.@, together with the experi- 
mental data by Ishida et al. jjJJ on underdoped Bi- 
2212 (T c = 79K). Using W(T), Eq.(|), with A=0.05 
and C=0.3, we find good agreement of our theoreti- 
cal results with the experimental data between T c and 
To w 130 K. Note that the external magnetic field applied 
in NMR experiments increases the density of free vortices 
by np(B) ~ B / cf> |Q. As a result, the phase coherence 
length (and in turn W) acquires a magnetic field depen- 
dence, Q 2 = 4n 2 K {T) (n F (T) + n F (B)). For T > T c , 



and the parameter set used above, we find that a mag- 
netic field of B — 10T increases the vortex density by 
nF(B)/riF{Tc) ~ 0.1; its contribution to W can thus 
be neglected. Thus the relaxation rate above T c should 
be independent of magnetic field for typical values of B, 
which is consistent with recent experiments by Gorny et 
al. @. 

In the above scenario, we neglected the effect of longi- 
tudinal phase fluctuations which arise from spin- wave like 
excitations. This is justified since their excitation spec- 
trum is very likely gapped by the Anderson-Higgs mech- 
anism and they are, consequently, irrelevant for the 
low-frequency thermodynamic properties of the under- 
doped cuprates. It was recently proposed in Ref. |22[ that 
longitudinal phase fluctuations are responsible for the lin- 
ear temperature dependence of the superfluid density at 
T « T c . FM pointed out that longitudinal phase fluctu- 
ations at T <C T c lead to a Wi onq ~ T. In this case it 



follows from Eq.(jFJ) that, for 63 Cu and 17 0, 1/T X T ~ T 
at T <C T c , in contrast to the experimentally observed 
1/T X T ~ T 2 §|. This result suggests that longitudi- 
nal phase fluctuations are absent in the superconducting 
state. 

We assumed above, following the argument applied to 
STS and ARPES experiments j^j , that transverse phase 
fluctuations are static on the time-scale of INS and NMR 
experiments which allowed us to neglect the quantum 
dynamical nature of the vortices. While this assump- 
tion likely holds for "fast" probes like INS, ARPES and 
STS where the quasi-particles are coupled to phase fluc- 
tuations for short times, it might be less justified for 
the much "slower" NMR experiments. In this light, the 
agreement of our theoretical NMR results with the exper- 
imental data, Fig. ||| is remarkable. However, the effects 
of the vortex quantum dynamics on various experimen- 
tal probes is still an open question which requires further 
study. 

In summary we propose a scenario for INS and NMR 
experiments in the pseudogap region of the underdoped 
cuprates. We argue that phase fluctuations of the su- 
perconducting order parameter drastically affect the fre- 
quency dependence of the spin susceptibility and can thus 
qualitatively account for the temperature dependence of 
the resonance peak. Moreover, we show that the spin- 
lattice relaxation rate, 1/T X T, is a direct probe for the 
strength of the phase fluctuations, as reflected in W(T). 
Finally, we showed that W(T) obtained from high fre- 
quency transport measurements is in good qualitative 
agreement with that extracted from STS experiments. 

It is our pleasure to thank A. H. Castro Neto, A.V. 
Chubukov, A. J. Leggett, A. J. Millis, D. Pines, R. Ra- 
mazashvili, J. Schmalian, R. Stern and M. Turlakov for 
valuable discussions. This work has been supported by 
Fundagdo de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP), the Center of Advanced Studies (CAS) of the 
University of Illinois (H.W.) and in part by the Science 
and Technology Center for Superconductivity through 
NSF-grant DMR91-20000 and by DOE at Los Alamos 



4 



(D.K.M). 



V.J. Emery and S.A. Kivelson, Nature 374, 434 (1995); 

ibid., Phys. Rev. Lett. 74, 3253 (1995). 

M. Franz and A.J. Millis, Phys. Rev. B, 58, 14572 (1998). 

H.-J. Kwon and A.T. Dorsey, Phys. Rev. B, 59, 6438 

(1999). 

A.V. Chubukov and J. Schmalian, Phys. Rev. B 57, 
R11085 (1998). 

V.G. Geshkenbein, L.B. Ioffe, and A.I. Larkin, Phys. 

Rev. B 55, 3173 (1997). 

S.C. Zhang, Science 275, 1089 (1997). 

P.A. Lee and X.-G. Wen, Phys. Rev. Lett. 78 4111 

(1997). 

H. Ding et al, Phys. Rev. Lett., 78, 2628 (1997). 
C. H. Renner et al, Phys. Rev. Lett, 80, 149 (1998). 
P. Dai et al, Science 284, 13 44 (1999); H.F. Fong et al, 



preprint, :ond-mat/9910041 



M. Takigawa, Phys. Rev. B 43, 247 (1991). 

K. Ishida et al, Phys. Rev. B 58, R5960 (1998). 

J. Corson, et al, Nature 398, 221 (1999). 

P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

D.K. Morr and D. Pines, Phys. Rev. Lett. 81, 1086 

(1998). 

P. Monthoux and D. Pines, Phys. Rev B 47, 6069 (1993). 
A.V. Chubukov, Europhys. Lett. 44, 655 (1997). 

D. K. Morr and R. Wortis, Phys. Rev. B 61, R882 (2000). 
We assume that A sw and c sw only exhibit a weak tem- 
perature dependence in the pseudo-gap region pof . 

P. Bourges in High Temperature Superconductivity, edited 
by S.E. Barnes et al. 207-212 (CP483 American Institute 
of Physics, Amsterdam, 1999) . 
K. Gorny et al, Phys. Rev. Lett. 82, 177 (1999). 

E. W. Carlson, et al, Phys. Rev. Lett, 83, 612 (1999). 
J. A. Martindale et al, Phys. Rev. B 47, 9155 (1993). 



5 



